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Abstract 

The temperature dependencies of specific heat and spin susceptibihty of a 
coupled d^2_y2 +idxy superconductor in the presence of a weak dxy component 
are investigated in the tight-binding model (1) on square lattice and (2) on a 
lattice with orthorhombic distortion. As the temperature is lowered past the 
critical temperature Tc, first a less ordered d^^-y^ superconductor is created, 
which changes to a more ordered dx2_y2 + idxy superconductor at Tci(< Tc). 
This manifests in two second order phase transitions identified by two jumps 
in specific heat at Tc and Td. The temperature dependencies of the supercon- 
ducting observables exhibit a change from power-law to exponential behavior 
as temperature is lowered below Td and confirm the new phase transition. 
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The unconventional high-Tc superconductors |T]] with a high critical temperature Tc have 
a complicated lattice structure with extended and/or mixed symmetry for the order pa- 
rameter . For many of these high-T^ materials, the order parameter exhibits anisotropic 
behavior. The detailed nature of anisotropy was thought to be typical to that of an extended 
s-wave, a pure ci-wave, or a mixed (s + exp(z6')(i)-wave type. Some high-Tc materials have 
singlet c?-wave Cooper pairs and the order parameter has d^2_y2 symmetry in two dimen- 
sions which has been supported by recent studies of temperature dependence of some 
superconducting observables 0-0]- In some cases there is the signature of an extended s- 
or (i-wave symmetry. The possibility of a mixed angular-momentum-state symmetry was 



suggested sometime ago by Ruckenstein et al. and Kotliar |jT^,|T5[. There are experimen- 



tal evidences based on Josephson supercurrent for tunneling between Pb and YBa2Cu307 



(YBCO) |jT6[|, and on photoemission studies on Bi2Sr2CaCu208+a; fl^l among others which 



are difficult to reconcile employing a pure s- or d-wave order parameter. These observations 



suggest that a mixed [s + exp(26')(i] symmetry is applicable in these cases |lT8| , p!9 



More recently Krishana et al. ||2^ reported a phase transition in the high-Tc supercon- 



ductor Bi2Sr2CaCu208 induced by a magnetic field from a study of the thermal conductivity 
as a function of temperature and magnetic field. Laughlin |^ have suggested that the new 



superconducting phase is the time-reversal symmetry breaking d^^^y^ + id^y state. Similar 



conclusion has been reached by Ghosh in a recent work from a study of the temperature 
dependence of thermal conductivity. This has led to the possibility of the transition to a 
(ia;2_y2 + id^y phase from a pure d^2_y2 phase of Bi2Sr2CaCu208. From a study of vortex 
in a (i-wave superconductor using a self-consistent Bogoliubov-de Gennes formalism, Franz 
and Tesanovic also predicted the possibility of the creation of a d.j.2_y2 + idxy super- 
conducting state. Although, the creation of the mixed superconducting state dx2_y2 + id^y 
is speculative, Franz and Tesanovic conclude that a dramatic change should be observed 
(in the observables of the superconductor) as the superconductor undergoes a phase tran- 
sition from a d.j.2^y2 state to a dx2-y2 + id^y state. In this work we study the effect of this 



phase transition on superconducting specific heat and spin susceptibihty in the absence of 
magnetic field. The general trend of observables under the (i^2_j^2- to dx'i-y2 + idr^y-wave 
phase transition of the superconductor, as studied here, is expected to be independent of 
the external magnetic field. We recall that there have been several studies on the formation 



of a mixed (s + i(i)-wave superconducting state from a pure (i-wave state ||2^-p9|. 

First we study the temperature dependence of the order parameter of the mixed d^-^-y-i + 
idxy state within the Bardeen-Cooper-Schrieffer (BCS) model [^O|,0. The BCS model for 
a mixed dx2_y2 + id^y state becomes a coupled set of equations. The ratio of the strengths 
of the dx2_y2- and d^y-wave interactions should lie in a narrow region in order to have a 
coexisting dx2-y2- and dxy-wave phases in the case of dx2-y2 + idxy symmetry. As the dx2-y2- 
wave (d^^y-wave) interaction becomes stronger, the d^^^-wave ((ia.2_y2-wave) component of the 
order parameter quickly reduces and disappears and a pure (ia;2_j^2-wave (d^^y-wave) state 
emerges. 

The order parameter of each of dx2-y2 and d^y states has nodes on the Fermi surface and 
may change sign along the Fermi surface. The s-wave order parameter does not have this 
property. Because of this, many ci-wave superconducting observables have power-law depen- 
dence on temperature, whereas the s-wave observables exhibit exponential dependence. We 
find that in the present coupled dx2_y2+idxy state the order parameter does not exhibit nodes 
and change of sign along the Fermi surface and exhibits a typical s-wave like behavior. Con- 
sequently, the observables in the coupled dx2-y2 +idxy state does not exhibit typical rf-wave 
power-law dependence on temperature, but rather a typical s-wave exponential dependence. 

For a weaker dxy-^sve admixture, in the present study we establish in the two- 
dimensional tight-binding model (1) on square lattice and (2) on a lattice with orthorhombic 
distortion another second-order phase transition at T = T^i < Tc, where the superconduct- 
ing phase changes from a pure dx2-y2-waye state for T > Td to a mixed dx2-y2 + idxy-wave 
state for T < T^. The specific heat exhibits two jumps at the transition points T = T^i 
and T = Tc- The temperature dependencies of the superconducting specific heat and spin 
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susceptibility change drastically at T = Tci from power-law behavior for T > to expo- 
nential behavior for T < Td. We find that the observables for the normal state are closer to 
those for a pure superconducting d^y state than to those for a pure superconducting (ij;2_j^2 
state. Consequently, superconductivity in dx2-y2 wave is more pronounced than in pure 
dxy wave. Hence as temperature decreases the system passes from the normal state to a 
"less" superconducting d^2_j^2-wave state at T — Tc and then to a "more" superconducting 
state dx2-y2 + idxy with dominating s-wave behavior at T = Tci signaling a second phase 
transition. 

The profound change in the nature of the superconducting state at T = T^i becomes 
apparent from a study of the entropy. At a particular temperature the entropy for the 
normal state is larger than that for all superconducting states signaling an increase in order 
in the superconducting state. In the case of the present dx2-.y2 -\- id^y state we find that as 
the temperature is lowered past Td, the entropy of the superconducting d^^-yi + idxy state 
decreases very rapidly (not shown explicitly in this work) indicating the appearance of a 
more ordered superconducting phase and a second phase transition. 

We base the present study on the two-dimensional tight binding model which we describe 
below. This model is sufficiently general for considering mixed angular momentum states, 
with or without orthorhombic distortion, employing nearest and second-nearest-neighbor 
hopping integrals. The effective interaction in this case can be written as 



Here Vi and V2 are the couplings of effective d^^-y^- and dxy-^awe interactions, respectively. 
As we shall consider Cooper pairing and subsequent BCS condensation in both these waves 
the constants Vi and V2 will be taken to be positive corresponding to attractive interactions. 
In this case the quasiparticle dispersion relation is given by 



Vkq = — Vi (cos kx — P COS ky) (cos Qx — P cos Qy) 



V2(sin kx sin ky) (sin Qx sin qy). 



(1) 



ek = — 2i [cos kx + (3 cos ky — ^ cos kx cos ky] , 



(2) 



4 



where t and (3t are the nearest-neighbor hopping integrals along the in-plane a and b axes, 
respectively, and 7t/2 is the second-nearest-neighbor hopping integral. 

We consider the weak-coupling BCS model in two dimensions with dx2_y2+idxy symmetry. 
At a finite T, one has the following BCS equation 

^--EV'.^^tauh^ (3) 

with iJq = [(eq — EpY + |Aqp]^/^, where Ep is the Fermi energy and the Boltzmann 
constant. The order parameter Aq has the following anisotropic form: Aq = Ai(cosg2,' — 
P cos qy) + iA2 sin sin qy. Using the above form of Aq and potential (|l|), Eq. becomes 
the following coupled set of BCS equations 

V 2kBT ^ ^ 

A, = V,E^^^^^^^^^^^tanh^ (5) 
V 2Eq 2kBT ^ ' 

where the coupling is introduced through i?q. In Eqs. (1) and (|^) both the interactions Vi 

and V2 are assumed to be energy-independent constants for |eq — Ep\ < ksTp, and zero for 

|eq — Ep\ > ksTo, where ksTo is the usual Debye cutoff. 



The specific heat is given by |31 



where /q = 1/(1 + exp^Ef^/ksT)). The spin susceptibility x is defined by ||T^,[T3 



x(T) = %^EA(i-A) (7) 

q 

where /itv is the nuclear magneton. 

We solved the coupled set of equations (|^) and (j^) numerically by the method of iteration 
and calculated the gaps Ai and A2 at various temperatures for T < Tc. We have performed 
calculations (1) on a perfect square lattice and (2) in the presence of an orthorhombic dis- 
tortion with Debye cut off /c^Td = 0.02586 eV {To = 300 K) in both cases. The parameters 



for these two cases are the following: (1) Square lattice — (a) t = 0.2586 eV, /3 = 1, 7 = 0, 
Vi = 0.73t, and V2 = 6.8t, = 71 K, = 28 K; (b) t = 0.2586 eV, /5 = 1, 7 = 0, 
Vi = 0.73t, and V2 = 7.9t, = 71 K, T^i = 47 K; (2) Orthorhombic distortion - (a) 
t = 0.2586 eV, p = 0.95, and 7 = 0, 14 = 0.97t, and V2 = 6.5t, = 70 K, = 25 K; (b) 
t = 0.2586 eV, p = 0.95, and 7 = 0, K = 0.97t, and l^a = 8.0t, = 70 K, = 52 K. 
For a very weak (i2,2„j^2-wave [d^y-wave) coupling the only possible solution corresponds to 
A2 = (Ai = 0). 

In Figs. 1 and 2 we plot the temperature dependencies of different A's for the following 
two sets of dx2_y2 +idxy-v^aiVe corresponding to models 1 and 2 above (full line — models 1(a) 
and 2(a); dashed line — models 1(b) and 2(b)), respectively. In both cases the temperature 
dependence of the A's are very similar. In the coupled dx2-y2 + id^y-wave as temperature 
is lowered past T^, the parameter Ai increases up to T = T^. With further reduction of 
temperature, the parameter A2 becomes nonzero and begins to increase and eventually both 
Ai and A2 first increases and then saturates as temperature tends to zero. Recently, the 
temperature dependencies of the order parameter of the d^^-y^ + is-wave superconducting 
state has been studied, where at T = T^, the transition from d^^-yi to dx2-y2+is state takes 
place. In that case, below T = T^i the dx'2-y2-wme component of the order parameter is 
suppressed, as the s-wave component becomes nonzero. No such suppression of the dx2_y2- 
wave takes place in this the dxy component appears. 

Now we study the temperature dependence of specific heat in some detail. The different 
superconducting and normal specific heats are plotted in Figs. 3 and 4 for square lattice 
[models 1(a) and 1(b)] and orthorhombic distortion [models 2(a) and 2(b)], respectively. In 
both cases the specific heat exhibits two jumps — one at Tc and another at T^. From Eq. 
(|) and Figs. 1 and 2 we see that the temperature derivative of |Aqp has discontinuities at 
Tc and Td due to the vanishing of Ai and A2, respectively, responsible for the two jumps 
in specific heat. For a pure dx2^y2 wave we find that the specific heat exhibits a power- 
law dependence on temperature. However, the exponent of this dependence varies with 
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temperature. For small T the exponent is approximately 2.5, and for large T {T ^ Tc) it 
is nearly 2. In the dx2-y2 + id^y model, for Tc > T > the specific heat exhibits dx2-y2- 
wave power-law behavior; for T < T^i the specific heat exhibits an s-wave like exponential 
behavior. For the d-wave model d^^-y^, Cs{Tc)/Cn{Tc) is a function of Tc and jS. In Figs. 3 
and 4 this ratio, for Tc = 70 K, is approximately 3 (2.5) for (3=1 (0.95). In a continuum 
(i-wave calculation this ratio was 2 in the absence of a van Hove singularity P^JTB| . We 



also calculated the specific heat for the pure d^y case. For square lattice with Vi = 9.0 we 
obtain = 67 K and C,(T,)/C„(rj = 1.82 and C,(T)/C„(r,) ~ {T/Tcf-^ for the whole 
temperature range. For orthorhombic distortion (3 = 0.95, with Vi = 9.0 we obtain Tc = 69 
K and Cs{Tc)/Cn{Tc) = 1.94 and C,(T)/C„(T,) ~ {T/Tc)^-^ for the whole temperature 
range. This power-law behavior with temperature in both the d waves is destroyed in the 
coupled dx2_y2 + idxy wave and for T < T^i, we find an s- wave-like exponential behavior in 
both cases. In both the uncoupled d waves the order parameter A has nodes on the Fermi 
surface and changes sign and this property is destroyed in the coupled d^^-yi + id^y wave, 
where the order parameter has a typical s wave behavior. 

In Fig. 5 we study the jump AC in the specific heat at Tc for pure s- and d-wave 
superconductors as a function of Tc, where we plot the ratio AC / Cn(Tc) versus Tc. For a BCS 
superconductor in the continuum AC/Cn{Tc) = 1.43 (1.0) for s-wave (rf-wave) independent 



of Tc [|T^,|T^,^ . Because of the presence of the van Hove singularity in the present model this 
ratio increases with Tc as can be seen in Fig. 5. For a fixed Tc, the ratio AC/Cn{Tc) is larger 
for square lattice (/? = 1) than that for a lattice with orthorhombic distortion {(3 = 0.95) for 
both s and d^^-y^ waves. However, for a ci^^y-wave superconductor AC/Cn{Tc) is smaller for 
square lattice (/5 = 1) than that for a lattice with orthorhombic distortion (/3 = 0.95). The 
jump in d^y wave is smaller than that for s and d^^^y^ waves. At Tc = 100 K, in the s-wave 
((ia;2_j^2-wave) square lattice case this ratio could be as high as 3.63 (2.92), whereas for dxy 
wave this ratio at 100 K is 1.15 (1.25) for square lattice (orthorhombic distortion). 

Next we study the temperature dependencies of spin susceptibility for square lattice and 



in the presence of orthorhombic distortion which we exhibit in Figs. 6 and 7, respectively. 
There we plot the results for pure dx^-y^, d^y, and s waves for comparison, in addition to 
those for models 1(a), 1(b), 2(a) and 2(b). In all cases reported in these figures ~ 70 
K. For pure ci-wave cases we obtain power-law dependencies on temperature. The exponent 
for this power-law scaling is independent of critical temperature T^. but vary from a square 
lattice to that with an orthorhombic distortion. For dx2_y2 wave, the exponent for square 
lattice (orthorhombic distortion, f3 = 0.95) is 2.6 (2.4). For dxy wave, the exponent for 
square lattice (orthorhombic distortion, (3 = 0.95) is 1.1 (1.6). For the mixed dx2-y2 + id^y 
wave, dx2^y2-wa:ve power-law behavior is obtained for Tc > T > T^. For T < T^, one has a 
typical s-wave behavior. 

In conclusion, we have studied the {d^i-yi + idxy)-^8cve superconductivity employing 
a two-dimensional tight binding BCS model on square lattice and also for orthorhombic 
distortion. We have kept the potential couplings in such a domain that a coupled {dx2-y2 + 
idxyY^SNe solution is allowed. For a weaker d^y admixture, as temperature is lowered 
past the first critical temperature Tc, a weaker (less ordered) superconducting phase is 
created in dx2-y2 wave, which changes to a stronger (more ordered) superconducting phase 
in {dx2-y2 + idxy) wave at T^. The dx2-y2 + idxy-wave state is similar to an s-wave-type 
state with no node in the order parameter. The phase transition at from a dx2_y2 wave 
to a dx2-y2 + idxy wave is marked by power-law (exponential) temperature dependencies of 
specific heat and spin susceptibility for T > Td (< Td). Similar behavior has been observed 
for a dx2^y2 + is-wave state [|2^. As the mixed state is s-wave like in both cases, from 
the present study it would not be possible to identify the proper symmetry of the order 
parameter — dx2-y2 + idxy opposed to dx2-y2 + is — and further phase sensitive tests of 
pairing symmetry in cuprate superconductors is needed. 

We thank Conselho Nacional de Desenvolvimento Cientifico e Tecnologico and Fundagao 
de Amparo a Pesquisa do Estado de Sao Paulo for financial support. 
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Figure Captions: 

1. The order parameters Ai, A2 in Kelvin at different temperatures for dx2_y2+idxy-waye 
models 1(a) (full line) and 1(b) (dashed Hue) for square lattice. 

2. The order parameters Ai, A2 in Kelvin at different temperatures for dx^-y^ + idxy- 
wave models 2(a) (full line) and 2(b) (dashed line) in presence of orthorhombic distortion 
(/? = 0.95). 

3. Specific heat ratio C{T)/Cn{Tc) versus T/Tc for models 1(a) and 1(b) for square 
lattice: 1(a) (full hue), 1(b) (dashed hue), d^y (dotted line), normal (dashed-dotted line). 
In all cases Tc ~ 70 K. 

4. Specific heat ratio C{T)/Cn{Tc) versus T/Tc for models 2(a) and 2(b) for orthorhombic 
distortion: 2(a) (full hne), 2(b) (dashed hue), d^y (dotted line), normal (dashed-dotted line). 
In all cases ~ 70 K. 

5. Specific heat jump for different Tc for pure s and d waves: s wave (solid line, square 
lattice), s wave (dashed line, orthorhombic distortion), d^^-yi wave (dashed-dotted line, 
square lattice), d^^-y^ wave (dashed-double-dotted line, orthorhombic distortion), d^y wave 
(dotted line, square lattice), dxy wave (dashed-triple-dotted line, orthorhombic distortion). 

6. Susceptibility ratio xiT)/x{Tc) for square lattice versus T/Tc- pure d^-^-y-i wave (solid 
line), pure d^y wave (dashed line), pure s wave (dashed-dotted line), model 1(a) (dotted 
line), model 1(b) (dashed-double-dotted line). In all cases Tc ~ 70 K. 

7. Susceptibility ratio x{T)/x{Tc) in presence of orthorhombic distortion versus T/Tc- 
pure dx2^y2 wave (solid line), pure d^y wave (dashed line), pure s wave (dashed-dotted line), 
model 2(a) (dotted line), model 2(b) (dashed-double-dotted hne). In all cases Tc ^ 70 K. 
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Figure 4 
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Figure 7 
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